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Abstract 

We study numerically the nature of the diffusion process on a hon- 
eycomb and a quasi-lattice, where a point particle, moving along the 
bonds of the lattice, scatters from randomly placed scatterers on the 
lattice sites according to strictly deterministic rules. For the honey- 
comb lattice fully occupied by fixed rotators two (symmetric) isolated 
critical points appear to be present, with the same hyperscaling re- 
lation as for the square and the triangular lattices. No such points 
appear to exist for the quasi-lattice. A comprehensive comparison is 
made with the behavior on the previously studied square and trian- 
gular lattices. A great variety of diffusive behavior is found ranging 
from propagation, super-diffusion, normal, quasi-normal, anomalous 
to absence of diffusion. The infiuence of the scattering rules as well 
as of the lattice structure on the diffusive behavior of a point particle 
moving on the all lattices studied so far, is summarized. 

KEY WORDS: Diffusion, Lorentz lattice gas, cellular automata, hon- 
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1 Introduction 

In a number of previous papers we have studied numerically the diffusion in Lorentz 
Lattice Gas Cellular Automata (LLGCA) for a variety of strictly deterministic scatter- 
ing rules^^"''] (for reviews see refs.[I], [4] and [6]). In particular new types of diffusion 
were found on the square and triangular latticest^^. In order to see to what extent the 
nature of the diffusion process depends on the type of lattice, we report here results 
for the honeycomb and quasi-lattices. Although some preliminary investigations were 
also made on the random lattice, they will not be considered here since they were not 
elaborate enough to be used for comparisons^]. The random lattice will be discussed 
in a later publication^*]. Since not all lattices considered in this paper are regular 
lattices, the results obtained allow us also a comparison of the diffusive behavior on 
a lattice in its dependence on the lattice structure as well as on the scattering rules. 

In LLGCA a point particle moves along the bonds of a lattice, whose sites are 
randomly occupied by stationary (in position) scatterers, which scatter the particle 
according to deterministic rules. We have studied two models for the scatterers, with 
about the simplest nontrivial scattering rules one can think of: either with right and 
left mirrors or with right and left rotators and considered two cases for each model: 
the scatterers do or do not change in character at collision, respectively. In the 
first case we consider fixed scatterers. i.e., stationary both in position and character 
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(discussed in ref. [1]), in the second case we consider flipping scatterers, i.e., stationary 
in position, but changing from right to left and vice-versa after a collision. 

In the case of fixed scatterers one can consider one particle or many particles 
moving simultaneously but independently of each other through the scatterers. For 
hipping scatterers, however, the motion of more than one particle through the scat- 
terers differs qualitatively from that of one single particle, since the hipping of the 
scatterers introduces (indirect) interactions between the particles as well as the scat- 
terers. We confine ourselves here to the motion of a single particle. 

Since straight lines do not pass through any lattice site on the honeycomb lattice 
and pass through some but not all lattice sites on the quasi-lattice, only the case 
of a honeycomb lattice and a quasi-lattice fully occupied by scatterers, i.e., with 
concentration of scatterers C = 1, has been considered, since otherwise a mixture of 
two models with different scattering rules would have to be introduced. 

The results of the investigations reported in this paper, which were carried out to 
over a million or more time steps, can be summarized as follows. 

1. Like on the square (for the fixed rotator model) and the triangular (for both the 
fixed rotator and the fixed mirror models) lattices, there appear to exist critical points 
for the fixed rotator model on the honeycomb lattice at which the diffusive behavior is 
anomalous and the trajectories obey a hyperscaling relation, while there is an absence 
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of difFusion-to which we will refer as no-diffusion-for all other concentrations. There 
is no-diffusion for any concentration of left or right fixed rotators on the quasi-lattice. 

2. There is no-diffusion for hipping scatterers (mirrors or rotators) on both the 
honeycomb and the quasi-lattices. 

3. A more detailed summary and comparison of the results for all lattices studied 
so far is given in Section 4. 

The organization of this paper is as follows. In Section 2, we study the honeycomb 
lattice for the two models mentioned above, by specifying the equations of motion 
for the moving particle. The Boltzmann approximation to the diffusion coefficient is 
given, the computer simulation method is described, and the results of the simula- 
tions for both fixed and fiipping scatterers are discussed. Section 3 discusses how to 
construct a Fibonacci quasi-lattice as well as the equations of motion for a particle 
moving on this lattice. The Boltzmann approximation is given, and the computer 
simulations as well as their results are discussed. In Section 4 we give a summary 
and a comparison with previous work. 

2 Honeycomb Lattice 

A. Equations of Motion 

The honeycomb lattice, shown in fig. I, has two kinds of sites which, after a rotation 
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over TT, transform into each other (cf.fig.2a) and play an identical role in the diffusion 
process. The scattering rules for the mirror and the rotator models are shown in 
fig.2b-c^. They lead to the following equations of motion for a particle for both fixed 
and fiipping rotators: 

ni{f+ei,t + l) = n2{r,t)mR{r,t) + ne{r,t)mL{r,t) 

n,{r + e,,t + l) = n,+i{f,t)mR{f,t) + n,_i{f,t)mL{r,t) (z = 2,3,4,5) 
ne{f+ee,t + l) = ni{r,t)mR{r,t) + n^{r,t)mL{r,t) (2.f) 

while for both fixed and fiipping mirrors, one obtains: 

ni{r + ei,t + I) = n2{r,t)mR{r,t) + ne{r,t)mL{r,t) 

n,{r + e,,t + l) = n,_i{f,t)mR{f,t) + n,+i{f,t)mL{r,t) (z = 2,4) 

n,{f+e,,t + l) = n,_i{f,t)mL{r,t) + n,+i{f,t)mR{f,t) (z = 3,5) 

ne{f+ee,t + l) = ni{r,t)mL{r,t) + n^{r,t)mR{r,t) (2.2) 

Here, the are the unit vectors along the six velocity directions i = f , 6 (see fig. 2); 
ni[f^t) = f or (z = f,...,6), if a particle is or is not at the lattice site r at time 
t, respectively; niR^L^r^t) = f or 0, if a right (left) rotator (in the rotator model) or 
a right (left) mirror (in the mirror model) is or is not at the lattice site r at time t. 



^We remark that for the fully occupied triangular lattice, if the moving particle turns 
±7r/3 upon collision with a scatterer, the particle moves on the honeycomb lattice consid- 
ered herein] . 



respectively. niR^r^t) or niL^r^t) are independent of t for fixed scatterers, while they 
depend on t for flipping scatterers. 
B. Boltzmann Approximation 

The Boltzmann approximation to the eqs.(2.f) - (2.2) is obtained by averaging 
both sides of the equations over all possible random configurations of the scatterer and 
the particle and ignoring any correlations between scatterer occupation and particle 
velocity at each lattice site. This leads for all cases, i.e., for fixed and fiipping mirrors 
or rotators to an equation of the form: 

/,(t + f) = /,(t) + X:r.,/,(t) (z = f,...,6) (2.3) 

where fiit) is the probability to find a particle with the velocity direction along at 
time t and Ti^ are the elements of the collision matrix T, given by: 
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for the rotator case, where Cl,r =< rnL^R^r^t) >, the concentration of left (right) 

scatterers on the lattice, respectively. 

Using the formula of Ernst and Bindert^^, we have: 
1 



(2.4) 



with 



(2.5) 



and (yi\ the normalized 1-component of the velocity of the particle in the basis: 

(2.6) 

For the matrix 2'™"'''°''^ the six eigenvalues are: 



/ | = ^n i -i -1 -i h 

y3^ '2' 2' ' 2' 2^ 



Ai = 0, A2 = -2, A3,5 = -1 + vr, A4,6 = -1 - vr 



(2.7) 



with six corresponding orthonormal eigenvectors: 
1 



(Ai 
(A2 
(A3 
(A4 
(As 
(Ae 



^6 



(1,1,1,1,1,1) 



-^(-1,1,-1,1,-1,1) 

1 2Cr-Cl ^ 2Cl - Cr _1 1 1, 

^/r 2^/V ' ' 2^/V ' 2' 2^/v' 2^ 

1 , Cl-2Cr ^ Cr - 2Cl _1 1 1. 

^ ' 2' 2^/v' 2' 



^3 



2yf 

1 Cr — Cl 1 - 

2Vr'"' 2Vr' 2' 2yf '2- 

Cl „ C*/? 1 C*/? — Cl 1 . 



^2Vf' ' 2yf'2' 2yf ' 2' 



(2.8) 
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where F = C£ — ClCr + C^. Using the projection operator ^ |Ai)(Ai| in (2.5), one 
obtains 



1=1 



imirror ^ 



while, similarly, for rotators one finds: 



(2.9) 



[rotator 



2{Cl + CI) 



(2.10) 



For the case of Cl = Cr = — , one has then 

Db = 1 (2.11) 

for both fixed and fiipping mirrors as well as rotators. 
C. Computer Simulations 

All the computer simulations were carried out on SiliconGraphics Indigo (SGI) 
with 32MB memory, SUN Sparc IPC with 24MB memory and VAX 3100 with 32MB 
memory. A virtual lattice of 90,000 X 90,000 sites was used. This lattice was con- 
ceptually divided into 300 X 300 blocks of 300 x 300 sites each. Two different arrays 
were used: one to record the position of a block, the other to record the position 
of the particle in the block. Actual memory was not assigned to a block (i.e., using 
"malloc" in C language) until the particle entered it. A fiag is used to mark whether 
a block has been visited by the moving particle or not: if the block was not visited, 
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i.e., the flag is 0, we assign memory to and put scatterers randomly on this block; 
otherwise, i.e., if the flag is 1, we continue to use the old conflguration of scatterers 
that was on the block before. After the particle flnished moving at the cut-off time 
step, we clear the memory and reset the flag for each block to 0. 

The advantage of this scheme^^^ is that only a small fraction of the memory for a 
90,000 X 90,000 array is actually used, since we do not have to reserve memory for 
those that are never visited by a particle. 

Another advantage of this memory allocation procedure appears when resetting 
the lattice to the blank condition, which must be done after each trajectory is com- 
pleted. For a large array this operation itself would take a signiflcant amount of time 
and would be required for each particle, even for very small closed orbits. In the 
present method, only the blocks that have been entered by the particle need to be 
reset when the trajectory is completed (i.e., using "free" in C language), signiflcantly 
reducing the average time required for this operation. 

About 30,000 independent particles, initially placed randomly on the lattice, were 
studied. In the case of flipping scatterers one particle was studied at a time. The 
calculations were done up to 2^° to 2^^ time steps. The statistical errors were de- 
termined by doing the calculations in two steps: flrst an average was made over all 
10,000 particles, with a different random conflguration of the scatterers for each par- 
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tide, then further averages were computed over typically three runs, involving three 
samples of 10,000 particles for each. The standard deviations of the mean are plotted 
as the error bars of the data in the figures. If the error bar does not appear, the error 
bar is inside the symbol. 
C.l Fixed Scatterers 

Like on the square and triangular lattices the diffusion is non-Gaussian since there 
are closed orbits on the honeycomb lattice for both the mirror and the rotator models. 

a. Cl = Cr = - 

In this case the two models behave identically, since they can be mapped into 
each other in a similar way as on the square latticet^^^. There is no-diffusion (Class 
IV), i.e., for increasing t the diffusion coefficient D[t) goes to zero, see fig. 3a, means 
that the mean square displacement A(t) is bounded, i.e., there are no extended closed 
orbits; consequently, the distribution function P(r, t) does not correspond to that of a 
Gaussian diffusion process, but exhibits a sharp peak near the origin at r^ax — 1 — 2, 
and its shape does not appear to change any more after about 2^^ time steps (cf.fig.3b). 
Since the behavior of the two models is the same, we show that for the rotator model 
only. 

b. Cl^Cr{Cl + Cr = 1) 

Like on the fully occupied square lattice, the rotator model on the honeycomb 
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lattice behaves very differently Irom the mirror model. For, while the mirror model 
exhibits no-diffusion (Class IV) at least for the concentrations 0.35 < Cl,r < 0.65, 
which we investigated (cf.fig.4a), the rotator model appears to possess two symmetric 
isolated critical points at CLcr,RcT — 0.541, Cr^^^Lct — 0.459, where the diffusive 
behavior is anomalous (Class II). This differs from the results reported by Catala et 
al.t^^l, who found a critical line from about Cl = 0.541 and C'r = 0.459 to C'l = 0.459 
and Cr = 0.541, of which only the end points are consistent with our results. For all 
other concentrations, away from these critical concentrations, the diffusive behavior 
is no-diffusion (Class IV) (see figs. 4b and 5). This behavior can also be seen in the 
difference in the number of open orbits as a function of time for Cl = Cl^^ 
Cl 7^ Clc^ (and similarly for Cr^^) (see fig. 6a). For the rotator model in fig. 6a, the 
number of open orbits No{t) has a maximum value at Cl^^ {^^ ^Rcr) approaches 
a behavior as t~^^'^^^\ i.e., we determined the size distribution exponent of open orbits 
characterized by r = y from No{t) ~ = t'^^"^ . At the same time, we determined 
independently the fractal dimension dj = ^ from the product of probability of an 
open orbit Po{t) and the mean square displacement of a particle on an open orbit 
Ao{t) divided by t, i.e., from Po{t)Ao{t)/t ~ r^/^t^/'^z/t which appears to approach a 
constant (cf.fig.6c). This suggests that the hyperscaling relation r — 1 = 2/df is valid 
at the critical points^^]. At all other concentrations, the orbits appear to close the 
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quicker the further away Cl (Cr) are from Cl^^ i^Rcr)- Although we expect that at 
Cicr (oi" CRct)} orbits will close eventually, it may take an infinite number of time 
steps to do so, i.e., there exist extended closed orbits^^]. We should point out that it 
appears to take much longer for the diffusion process on the honeycomb lattice at a 
critical point to reach its asymptotic behavior than on the square and especially the 
triangular lattices. To what extent this is related to the small coordination number 
(3) of the honeycomb lattice is unclear. We are therefore not as certain about the 
existence of the two critical points and the hyperscaling relation for the honeycomb 
lattice, as we are for the other two lattices. In order to settle this point unambiguously 
we would have to extend our calculations far beyond our present maximum of 2^^ time 
steps, where each data point in figs. 4a, 6a and 6c typically already takes a few weeks 
on our SGI. 

For the mirror model, no such critical points (cf.fig.6b) are found. In fact, the 
number of open orbits decreases gradually with increasing values of Cl when Cl 
varies from to 0.5, where it is a minimum and similarly for < Cr < 0.5. This is 
because it is more difficult in the mirror model for a moving particle to make a closed 
orbit when there are more mirrors of one type than the other, due to an increased 
possibility of zig-zag motion (cf.fig.2c in ref. [1]). 
C.2 Flipping Scatterers 
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Like on the fully occupied square and triangular lattices, the diffusive behavior of 

the mirror and the rotator models on the fully occupied honeycomb lattice is identical. 

This can be argued in the same way as for fixed scatterers at Cl = Cr = |- We find 

that the diffusion coefficient goes to zero and that the radial distribution function 

P(r, t) does not appear to change anymore after a number of time steps tcr ~ 2^^, 

see fig.7a-b (in fig. 7b, since the behavior of the two models is the same, we show 

that for the rotator model only), while the maximum of P(r, t) appears to remain 

fixed at a distance r = r^ax ~ 10 lattice distances from the origin. We classify this 

diffusive behavior still as Class IV, since it is similar to that found for fixed scatterers 

in Section C.l, except that r^ax is much larger now. We observe that after about 

2^^ time steps virtually all particles are in closed orbits (see fig. 8), suggesting that 

all trajectories will eventually become closed orbits. The smallest closed orbit can be 

considered to consist of two "refiectors" (i.e., the closed orbit consists of two parts 

connected by a single line, where each of these two parts is called a "refiector")^ with 

a period of 92 time steps, where the particle suffers 92 collisions on 34 lattice sites 

(cf.fig.9). This closed orbit is much larger than the corresponding one of 6 lattice 

sites for fixed scatterers, which is the origin of the larger value of r^ax for the hipping 
^L. Bunimovich and S. Troubetzkoy use a more general definition of a "reflector" I^^l. 
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scatterer model. 

b. Cl^Cr{Cl + Cr = 1) 

In this case, the diffusive behavior ol both the rotator and the mirror model is 
Class IV, but it is not entirely identical for the two models. 

For the flipping rotator model, the diffusion coeflicient takes a shorter time for 
Cl 7^ Cr to vanish than for Cl = Cr (see flg.7a). This is so, because for Cl 7^ Cr, 
when one type of rotators is more numerous than the other, the closed orbits are less 
extended than when Cl = Cr (cf.flgs.lOa-b), since for Cl 7^ Cr the tendency for the 
moving particle to bend in one direction and then return to its original position is 
larger than when an equal number of left and right rotators is present. 

We remark that the closed orbits found in our simulations for Cl 7^ Cr appear 
to be of a different type than those observed for Cl = Cr, viz. without "reflectors". 
This difference is surely only apparently so, since for all Cl and Cr 7^ closed orbits 
with or without "reflectors" will occur. However, their frequency of occurrence will 
depend very much on Cl/Cr. Thus, while for Cl = Cr those with "reflectors" are 
the only ones seen in the simulations, for Cl 7^ Cr other types of orbits appear as 
well (see flg.lOa). Similarly on the square lattice for flipping rotators, "reflectors" 
are not necessary for closed orbits either (cf.flgs.lOc-d). However, on the triangular 
lattice, for both flipping rotators and mirrors, closed orbits never contain "reflectors" 
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(as defined above) (see figs.lOe-f). This can be understood as follows: for a "reflector" 
to exist, it is necessary that the the moving particle goes through at least one bond 
at least once in opposite directions, i.e., the moving particle must turn over an angle 
of an odd number times of tt, consistent with the scattering rules. However, the latter 
is not possible for the triangular lattice because an angle of an odd number times of 
TT can not be produced by any combination of the possible scattering angles ib27r/3 
and 0. 

For the flipping mirror model, the diffusion coefficient goes to zero for both Cl 7^ 
Cr and Cl = Cr. However, for this model it takes a longer time for the diffusion 
coefficient to vanish when Cl 7^ Cr than for Cl = Cr^ see fig. 7a, since there exists 
now a tendency for the particle to propagate (zig-zag motion) when the particle hits 
a region with more right (left) mirrors than left (right) mirrors. 

3 Quasi-lattice 

A. Introduction 

The quasi-lattice we studied is a Fibonacci lattice rather than a Penrose lattice, 
since the former gives a higher density of vertices (for a given lattice size). It is 
shown in fig.ll. It was constructed in the following wa,y^^^\ First a "star" of 5 two- 
dimensional vectors Si, S2, ■Ss, ■5*4, ■55 is drawn (cf.fig.l2a). This "star" has pentagonal 
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orientational symmetry, i.e., the angle between each pair ol adjacent vectors is 27r/5. 
Next a grid, i.e., a set ol, in principle, infinite, but here in lact 73, quasi-periodically 
spaced parallel lines are introduced perpendicular to each star vector Si[i = 1, ...,5) 
with spacings: 

1 Tl 

Xn,z = rn,t ■ s, = T,[n + a,-\ [ h AJ] (3.1) 

Here Xn,i is the distance along the direction Si [i = 1,...,5) between the origin and 
the n-th parallel line; fn^i is a point on the n-th parallel line perpendicular to Si[i = 
1,...,5); [/J is the fioor function, denoting the integer part of / and n = 1,2,....,A^ 
(here N is 73). These 5 grids compose the grid-space. Furthermore, Ti^ai, pi^cTi and 
fii are constant parameters, of which the (Ti must be irrational real numbers in order to 
get non-periodic spacings"^. Eq.(3.1) defines a quasi-periodic "sequence of intervals" 
such that the intervals Ax^ = Xn^i — Xn-i,i between two adjacent parallel lines have 
the property: 

. _ / if In/a, + /3J - [(n - l)/a, + /3J = 0, 

' ~ 1 r,(l + 1/p,) if [n/a, + /3J - L(n - l)/a, + /3J = 1. ^''-'^ 

That is, there are only two possible intervals between adjacent parallel lines for 

each grid, Ti[l -\- l/pi) and T^, which appear in a quasi-periodic sequence, where 



^For, if we choose for ai in eq.(3.1) rational real numbers, e.g., ai = Pi/qi with both pi 
and Qi positive integers, we will get periodic spacings, Xm,i = T'i(m + ji) and therefore a 
periodic lattice with period T'i. Here m is an integer m = 1,2, [A^/kJ , where each m 
stands for a set of pi adjacent lines in the sequence of N lines, T'i = Ti(pi + Qi/pi) and 
7i = (tti + [I3i\/pi)/{pi + qi/pi)- 
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the ratio of the number of Ti[l + l/pi) distances to the number of Ti distances 
equals 1/(1 — (Ti). Thus, the parameter (Ti determines the relative frequencies of 
the two different spacings in the sequence and pi determines the ratio of the two 
spacing lengths that occur between the lines in each grid. In our calculations we 
took the parameters pi = (Ti = t = (a/5 + l)/2 {i = 1,...,5), the golden mean 
and we chose T, = 0.25, 0.8, 0.25, 0.7, 0.7; a, = -30.0, -30.0, -31.0, -72.0, -71.0 and 
l3i = 0.7,0.8,1.1,1.2,0.9, for i = 1,...,5, respectively, in order to make 5 grids that 
overlapped as much as possible. In this way a maximum number density of ver- 
tices of the quasi-lattice was obtained. The grid thus generated is called a Fibonacci 
pentagrid^^'*], to which we restrict ourselves here. The effect of a change in ai is simply 
to translate each entire grid in the direction Si over a distance Tj-ai, whereas a change 
in l3i alters the sequence of long and short spacings. In each grid, each line normal 
to the star vector Si is labeled by an integer ki which represents the location of its 
position along the Si direction. The lines divide the grid space into non-intersecting 
open regions through which no lines pass (the regions can be arbitrarily small). Each 
such region is specified (uniquely) by M (here 5) integers (A;i, A;2, kM)'- if Xq is any 
point in the region, then ki is the label of the line normal to Si such that Xq lies 
between the lines labeled by ki and ki -\- 1. Finally the generalized dual methodt^"^! is 
applied to find the dual of the pentagrid. The "dual" is constructed by mapping each 
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M 

open region in grid space into a point t = ''^kiSi, which Res in a two dimensional 

8 = 1 

space that we shall call the "cell-space". The points t are the vertices of a packing of 
the quasi-lattice by unit cells (rhombuses). So the dual transformation maps a grid- 
space into a cell-space such that open regions in grid-space are mapped into points in 
cell-space and the points in grid-space are mapped into the open regions (rhombuses) 
in cell-space. This gives a tiling of the plane with two differently shaped rhombuses 
(one is a thin rhombus with angles 36°, 144°, 36° and 144°, the other is a fat rhombus 
with angles 72°, 108°, 72° and 108°) and produces the quasi-lattice shown in fig. 11. 
The lengths of all bonds of the quasi-lattice are equal and are chosen to be unity. 
Since our quasi-lattice has no translational or rotational symmetries, we define two 
lattice sites as the same, if one lattice site can be mapped into the other by a pure 
translation. The Fibonacci quasi-lattice has 658 different lattice sites'*. If each lattice 
site were to appear on the lattice with the same probability, i.e., 1/658, then the 
average coordination number would be 5.6839. Since in our lattice, not every lattice 
site appears with the same probability the average coordination number is actually 
about 3.9887. 

B. Equations of Motion 

■*The lattice sites can differ in three ways: a) different number of bonds to the nearest 
neighbors; b) same number of bonds to the nearest neighbors, but different angles between 
adjacent bonds; c) same number of bonds to the nearest neighbors and same angles between 
adjacent bonds, but different orientation of these bonds. 
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We restrict ourselves to the rotator model, where the particle will turn to its left 
(right) over the largest available angle between lattice bonds if there is a left (right) 
rotator (cf. fig. 12c). This leads to the following equations of motion for both fixed and 
fiipping rotators: 

n,(f + e„t + 1) = [n,+i(f,t)r+iX^1 + n,+2{r,t)r+^\r) + 

[i = 1,...,10; mod 10) 

Here is the unit vector defining the velocity direction z (z = 1, 10) and ni(r, t), 
niR^r^t) and niL^r^t) have the same meaning as in the equations for the honeycomb 
lattice. P^[f) is a geometric factor, which can only have the values 1 or 0, depending 
on whether the particle can or cannot change from to e*,-, respectively, at r. In the 
square brackets before niR^r^t) or niL^r^t) in eq.(3.3), only one P^[f) is ^ 0, viz. 
that one which allows the particle coming in with velocity direction i to turn to the 
allowed direction j to its right or left, respectively. 
C. Boltzmann Approximation 

A Boltzmann approximation to the eq.(3.3) is obtained by averaging both sides of 
the equation over all possible random configurations of the scatterers and the particle 
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and ignoring any correlations between scatterer occupation and particle velocity at 
each lattice site. One obtains then the equations: 

10 

f^{t + i) = m + j2T,,m (z = i,...,io) (3.4) 

where fi{t) is the probability to find a particle with the velocity direction along at 
time t. Tij are the elements of a 10 x 10 collision matrix T, with off-diagonal matrix 
elements proportional to C, or 0% For C, = Cn = \, to which we restrict ourselves 
here, T reads: 



-2 < ^ > <I*'^> <I°^> </'-'> </''-'> </"-'> <I^°^> 



2 



\ 

<J-i^> _2 <I'^^> <I*'^> <I°^> < 1° ^ > <I''^> < I" ^ > < 1^°"^ > ^ 
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<I > <I > -2 <I*'^> <I°'^> <I°'^> <I''^> <I'''^> <I 

<jl''> <j21> <j34> -2 <f^*> <l7''> <l8''> <jlO' 

<jl5> <j2 5> <j3 5> <l''5> -2 <I^5> <I^5> <l3 5> <j9 5> q 

<l2''> <l3''> <r^'^> -2 <I^''> <l3''> <I^''> < ji"' 

< jl 7 > < l3 7 > < J* 7 > < l5 7 > < l6 7 > _2 < > < > < jlO 

< 3 > <j23> <l''3> <fi^> <fi^> < l'' ^ > -2 <I^3> <ll''i 

<jl9> <j29> <j39> <l59> <J«=9> <j79> <j89> -2 < jl" ' 

<jllO> <j2 10> <j3 10> <j4 10> <j6 10> <j7 10> <j8 10> <_r9 10> _2 



(3.5) 

Here < > is the fraction of all lattice sites for which r^[f) is not zero. The values of 
< > for the Fibonacci lattice are given in the Appendix. Using the method of Ernst 
and Bindert^l as above (cf.eqs.(2.4 - 2.5)) with the Tij of eq.(3.5), the diffusion coeffi- 
cient in the Boltzmann approximation can be obtained. Since the diffusion tensor is 

not isotropic, we use for the Boltzmann diffusion coefficient Db = -^{D^x^B + Dyy^s)- 

.1 TT 27r 2tt ^1 ^ 27r 27r TT, 

Usmg that {Vr\ = — ^( 1, cos — , cos — , — cos — , — cos — , —1, — cos — , — cos — , cos — , cos — 

^ ^ ' ^5^ 5' 5 ' 5 ' 5' ' 5' 5 ' 5 ' 5' 

1 27r TT TT 27r 27r tt tt 27r 

and {vy \ = — ^(0, cos — , cos — , cos — , cos — , 0, — cos — , — cos — , — cos — , — cos — ), 

^ ^' ^5^ 5 ' 5' 5' 5 ' ' 5 ' 5' 5' 5 ^' 
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which are the x and y components of the velocity of the particle in the [i = 1, 10) 
basis, respectively (cf.eq.(2.6)), leads then with eqs.(2.4) and (2.5) to a value Db — 
0.165. This value can be compared with the value Db — 0.1661 it 0.0005 as deter- 
mined by computer simulation for a probabilistic model with a scattering rule of equal 
probability for the particle to scatter over the largest angle to its left or to its right. 
D. Computer Simulations 

A unit cell of 22,350 lattice sites and Ziff et al.'st^l method described above for the 
honeycomb lattice were used. Since it is not possible here to impose strictly periodic 
boundary conditions, we chose the following quasi-periodic boundary conditions: each 
time that a particle leaves the unit cell along a bond with a given velocity at one side, 
it reappears on the opposite side with the same velocity direction along the bond 
closest to where it would have reappeared for periodic boundary conditions. About 
30,000 independent particles were used and the calculations were pursued up to 2^° 
time steps. For flipping scatterers one particle was studied at a time. As explained 
for the honeycomb lattice, the standard deviations of the mean are plotted as the 
error bars of the data in the figures. 
D.l Fixed Scatterers 

In this case, we found no-diffusion (Class IV) everywhere, since all orbits clearly 
seemed to be closed after a finite number of time steps. Therefore, unlike in the 
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honeycomb lattice, there appear to be no critical points. As a consequence, the 
mean square displacement is bounded, as can be seen this both from the diffusion 
coefficient and the number of closed orbits (see figs.f 3a-b). We also checked the effect 
of boundary conditions on the computer simulation results and found virtually no 
difference in the diffusive behavior between a basic block with 22,350 or f0,940 sites 
(cf.fig.f3c). 

D.2 Flipping Scatterers 

Like on the honeycomb lattice, the diffusion coefficient goes to zero, while the 
distribution function P(r, t) exhibits a sharp maximum at r = r^ax ~ fO lattice sites 
and does not appear to change anymore after 2^^ time steps (cf.figs.f 4a-b). It seems 
that all orbits will eventually be closed. Because of the irregularity of the lattice 
we have not been able to identify yet the smallest closed orbit; figs.f5a-b give two 
examples of closed orbits. 

4 Discussion 

f. We summarize the diffusive behavior for the fixed and fiipping scatterer models 

for Cl or Cr ^ on all four lattices investigated here as well as in ref. [f]^, in Table 

®For the diffusive behavior for C = Cl (or C = Cr), i.e., a lattice covered by only one 
type of scatterers, see ref.[6]. 
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I and II, respectively. 

The two tables reflect how the diffusion process differs on the fully occupied hon- 
eycomb and quasi-lattices from those studied before and enable us to get a sense of 
the influence of the lattice structure as well as of the scattering rules on the diffusion 
process of a particle on these lattices. We have the following comments. 

a) In Table I critical lines or points for the fixed rotator model occur for all 
lattices except the quasi-lattice. For the mirror model on the square lattice, the 
critical behavior occurs for all concentrations along the critical line C = 1. In Table 
II, a phase transition occurs at C = I for the rotator model on the square lattice as 
well as for both the rotator and mirror models on the triangular lattice. 

One can interpret the results in Tables I and II in terms of dynamical analogues 
of phase transitions, critical points and critical lines. A transition from Class IV 
to Class II exhibits some analogy to a second order phase transition, in that the 
behavior implies the appearance of extended closed orbits (cf. long range correlations 
at critical points or on critical lines) when before trapping, i.e., only short closed 
orbits, occurred (cf. short range correlations away from critical points or lines). On 
the other hand, the less subtle transitions from Class II to super-diffusion or Class 
IV to normal diffusion (Class I) or quasi-normal diffusion to propagation bear more 
resemblance to a first order phase transition. It is not clear whether this analogy can 
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be pushed beyond this descriptive stage. 

b) The triangular lattice is an exception in that the rotator and mirror models al- 
ways behave the same for any concentration of scatterers for both fixedt^^ and flipping 
scatterers for the same reason as argued in ref. [f] for fixed scatterers. 

c) In Tables I and II, the diffusive behavior of the moving particle covers a wide 
variety of different cases: propagation, super-diffusion, normal, quasi-normal, anoma- 
lous and no-diffusion. 

2. We note that the quasi-normal behavior for < C < f in the triangular lattice 
listed in Table II, differs from the normal behavior mentioned in ref. [2] before, because 
the finite probability for closed orbits to occur was overlooked in this paper^. This 
makes the diffusion process strictly speaking non-Gaussian[^°] . However, the fraction 
of closed orbits is so small (see fig. f 6a) that for all practical purposes the diffusion 
can be considered Gaussian; as a consequence we call it quasi-normal (cf.fig.f6b). We 
conjecture that a similar situation may obtain in d = 3 for the simple cubic lattice, 
for instance, where closed orbits will occur relatively frequently for short times when 
the particle hovers near the origin, but where, with increasing time, return to the 
origin, once the particle is away, will become increasingly difficult. 

3. We note that the mirror model for fixed scatterers is time reversal invariant 

®This was noticed by L. Bunimovich and S. Troubetzkoy (see ref.[10]) and independently 
by one of us (F. W., unpublished). 
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on the square and the triangular lattices but not on the honeycomb lattice. This is 
due to the small scattering angle (vr/3) of the particle on the honeycomb lattice as 
compared to those on the square and the triangular lattices, where they are 7r/2 and 
27r/3, respectively. However, the flipping mirror model on the honeycomb lattice is 
time reversal invariant, while this is not so on the other two lattices. 

4. As to the nature of closed orbits for flipping scatterers, empty sites are necessary 
on the square and triangular lattices, but not on the honeycomb and quasi-lattices 
(C = f). "Reflectors" occur but are not necessary for closed orbits on the square, 
honeycomb and quasi-lattices, i.e., closed orbits may or may not include "reflectors", 
while on the triangular lattice, there are no "reflectors" (cf.figs.f Oa-f ). 

5. Even when all orbits close, there are for all lattices a number of different cases 
to distinguish of how this closing occurs. 

(z). Many trajectories close very near the origin (where P has a maximum at 
r = ^max ~ f — 2) but most trajectories are extended and close gradually when 
t — )• oo, as shown in the slow approach of P[r^t) to its asymptotic shape (cf.figs.2a 
and 7a in ref. [f]). This obtains in the case of anomalous diffusion (Class II) for fixed 
scatterers. Then the decrease in the number of open orbits (characterized by r) and 
the mean square displacement (characterized by dj) are connected by the hyperscaling 
relation r — f = 2/df. 
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(a). Almost all orbits close very near the origin at r = r^ax ~ 1 — 2 and P[r^t) 
is virtually stationary for t > tcr (cf.fig.3b and fig. lib in ref. [f]). This obtains in 
the case of no-diffusion (Class IV) for fixed scatterers. All orbits close after a finite 
number of time steps. 

(Hi). Almost all orbits close near the origin at 10 and P(r,t) is 

virtually stationary for t > tcr ^ 2^^ (cf.figs.7b and 14b). This obtains in the case of 
no-diffusion (Class IV) for hipping scatterers. All orbits close after a finite number 
of time steps. 

(iv) . Only a small fraction of trajectories close near the origin, yet P[r^t) has a 
maximum at r^ax — 1 — 2 for all sufficiently large t. The particle appears to move 
overwhelmingly in unbounded trajectories (cf. fig. 16b). This obtains in the case of 
quasi-normal diffusion for both hipping rotators and mirrors on the triangular lattice 
for C < 1. 

(v) . All trajectories will close eventually with only a few closed orbits near the ori- 
gin, while the overwhelming majority carry particles whose mean square displacement 
grows faster than t (cf. fig. 7b in ref. [1]). This obtains in the case of super-diffusion 
for < C < 1 for the fixed mirror model on the square lattice. The probability to 
find a closed orbit decreases as ~ l/lntt^l. 

6. Although the theorems of Bunimovich and Troubetzkoy rigorously establish 
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some of the results found in these computer simulations on the square and the tri- 
angular lattices and we were able to present physical arguments for a number of 
the observed behaviors, clearly a deeper understanding of the diffusion process of 
a particle on a lattice in its dependence on lattice structure and scattering rules is 
lacking. 
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Figure Captions 



Fig.l. Part of a honeycomb lattice. 

Fig. 2. (a) The 6 velocity directions on the honeycomb lattice; after rotation over tt, 
lattice site I transforms into lattice site II and vice versa; (b) examples of rotator and 



Fig. 3. (a) Diffusion coefficient D as a function of the time t on a logio-logio scale for 
the fixed rotator model on the honeycomb lattice for Cl = Cr = 0.5; (b) correspond- 
ing radial distribution function P[r^t) as a function of distance r from the origin at 
t = (O), t = 2^1 (+) and t = 2^^ (□), respectively. 

Fig. 4. (a) Diffusion coefficient D as a function of the time t on a logio-logio scale 
for the fixed mirror model on the honeycomb lattice for Cl = 0.55, Cr = 0.45 (O), 
Cl = 0.6, Cr = 0.4 (+) and Cl = 0.65, Cr = 0.35 (□); (b) Diffusion coeflicient D 
as a function of the time t on a logio-logio scale for the fixed rotator model on the 
honeycomb lattice for Cl = Cr = 0.5 (O), Cl = 0.541, Cr = 0.459 (□) Cl = 0.55, 
Cr = 0.45 (+) and Cl = 0.6, Cr = 0.4 (x); critical concentration is printed in bold 



Fig. 5. Phase diagram for the fixed rotator model on the fully occupied honeycomb 
lattice. 




face. 
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Fig. 6. Number of open orbits out of 10,000 trajectories on the honeycomb lattice as 
a function of t on a log2-log2 scale for (a) fixed rotator model: Cl = Cr = 0.5 (O), 
C'l = 0.53, Cr = 0.47 (+), Cl = 0.541, Cr = 0.459 (□), C'l = 0.55, Cr = 0.45 
(x), Cl = 0.57, Cr = 0.43 (A) and Cl = 0.6, Cr = 0.4 (*); (b) fixed mirror model: 
Cl = Cr = 0.5 (O), Cl = 0.55, Cr = 0.45 (+), Cl = 0.6, Cr = 0.4 (□), Cl = 0.65, 
Cr = 0.35 (x); (c) contribution of diffusion coefficient from open orbits Po{t)Ao{t)/t 
as a function of time t on a log2-log2 scale for fixed rotator model on the honeycomb 
lattice for Cl = 0.541, Cr = 0.459 (O), Cl = 0.53, Cr = 0.47 (+) and Cl = 0.55, 
Cr = 0.45 (□). 

Fig. 7. (a) Diffusion coefficient D as a function of the time t on a logio-logio scale for 
the hipping rotator and mirror models on the honeycomb lattice for Cl = Cr = 0.5 
for hipping rotator (O) and fiipping mirror (+), Cl = 0.8 and Cr = 0.2 for fiipping 
rotator (□) and fiipping mirror (x); (b) corresponding P[r^t) as a function of r for 
the fiipping rotator model on the honeycomb lattice for Cl = Cr = 0.5 at t = 2^^ 
(O), t = 2^^ (+) and t = 2^^ (□), respectively. 

Fig. 8. Fraction of closed orbits on the honeycomb lattice for the fiipping rotator 
model (O) and fiipping mirror model (+) (indistinguishable) for Cl = Cr = 0.5, as 
a function of log2t. 

Fig. 9. An example of the smallest closed orbit for the fiipping rotator model on 
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the honeycomb lattice with period of 92 time steps on 34 lattice sites with two "re- 
flectors": A and B. Up-triangles stand for right rotators and down-triangles for left 
rotators, respectively. 

Fig. f 0. A few typical examples of closed orbits for the flipping rotator model on dif- 
ferent lattices. Up-triangles stand for right rotators; down-triangles for left rotators 
and squares for empty sites, respectively, (a) Period of 2209 time steps on 256 lattice 
sites without "reflectors" and empty sites for Cl = 0.95, Cr = 0.05 on the honeycomb 
lattice; (b) period of 852 time steps on f82 lattice sites with two "reflectors" (A and 
B) but without empty sites for Cl = Cr = 0.5 on the honeycomb lattice; it has a 
shorter period, but is more extended than that in fig. f 0(a) for Cl 7^ Cr; (c) period 
of 2697 time steps on 232 lattice sites without "refiectors" but with empty sites for 
Cl = 0.95, Cr = O.Of on the square lattice; (d) period of 337 time steps on 6f lattice 
sites with two "refiectors" (A and B) and empty sites for Cl = 0.95, Cr = O.Of on the 
square lattice; (e) the smallest closed orbit of period of f 8 time steps on 7 lattice sites 
without "refiectors" but with one empty site on the triangular lattice; (f) period of 
36 time steps on f3 lattice sites without "refiectors" but with empty sites for C < f .0 
on the triangular lattice. 
Fig. f f . Part of a Fibonacci quasi-lattice. 

Fig. f2. (a) Five- "star" with pentagonal orientation symmetry; (b) the fO velocity 
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directions on the Fibonacci quasi-lattice; (c) rotator scatterers. 

Fig. 13. (a) Diffusion coefficient D as a function of the time t on a iogio-iogio scaie 
for the fixed rotator modei on a Fibonacci quasi-iattice for Cl = Cr = 0.5 (O), 
Cl = 0.57, Cr = 0.43 (+) and Cl = 0.55, Cr = 0.45 (□); (b) as in (a) the number 
of closed orbits for Cl = Cr = 0.5 (O), Cl = 0.51, Cr = 0.49 (+), Cl = 0.52, 
Cr = 0.48 (□), Cl = 0.53, Cr = 0.47 (x), Cl = 0.54, Cr = 0.46 (A), Cl = 0.55, 
Cr = 0.45 (*), Cl = 0.56, Cr = 0.44 (O), Cl = 0.57, Cr = 0.43 (+) and Cl = 0.58, 
Cr = 0.42 (□); (c) as in (a) for Cl = Cr = 0.5 on a basic block with 22,350 sites (O) 
and a basic block with 10,940 sites (+) (indistinguishable). 

Fig. 14. (a) Diffusion coefficient D as a function of the time t on a logio-logio scale 
for the hipping rotator model on a Fibonacci quasi-lattice for Cl = Cr = 0.5 (O), 
Cl = 0.8, Cr = 0.2 (+) and Cl = 1.0, Cr = (□); (b) corresponding P{r,t) as 
a function of r for Cl = Cr = 0.5 at t = 2^^ (O), t = 2^^ (+) and t = 2^^ (□), 
respectively. 

Fig. 15. (a) An example of a closed orbits for the hipping rotator model on the 
Fibonacci quasi-lattice with period of 1129 time steps on 274 lattice sites with two 
"refiectors": A and B; (b) an example of a closed orbit with period of 5997 time steps 
on 348 lattice sites without "refiectors". Up-triangles stand for the right rotators and 
down-triangles for the left rotators. 
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Fig. 16. (a) Fraction of closed orbits for the flipping rotator model on the triangular 
lattice as a function of log2t for different concentrations of rotators for Cl = Cr; 
from top to bottom: C = 0.85, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1 (the last two are 
too small to see), respectively; (b) P(r, t) as a function of r for the flipping rotator 
model (circles) and flipping mirror model (diamonds) on the triangular lattice for 
C = 0.85 and Cl = Cr at t = 2^^. Smooth curves (indistinguishable) represent the 
corresponding Gaussian distribution (note the few points representing closed orbits 
near the origin). 
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Table Captions 



Table I: Fixed Scatterers. 
Table II: Flipping Scatterers. 



TABLE I - FIXED SCATTERERS 



1 , ?i "i" "i" 1 r*f^ 




\/1 1 T*T*r^T* 




except on two critical linest^^ — > Class II 


II ( 1 V '^11"I~\^T* /H 1 TT1 1 1 *~in 

U ^ ^ 1 ? O Lipci -Qlli LiolOn 


C = 1 ^ Class 11. 


triangular 


< C < 1 ^ Class IV 
except on one critical line 
Cl = Cj?[i] — > Class II 


< C < 1 ^ Class IV 
except on one critical line 
Cl = Cj?[i] — > Class II 


honeycomb 


C = 1 ^ Class IV 
except for two critical points > Class II 


C = 1 ^ Class IV 


quasi 


C = 1 ^ Class IV 
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TABLE II - FLIPPING SCATTERERS 



Lattice 


Rotator 


Mirror 


square 


< C < 1 ^ Class IV 


< C < 1 ^ Normal (Class I) 


C = 1 ^ Normal (Class I) 


triangular 


< C < 1 — > Quasi-normal 


< C < 1 — > Quasi-normal 


C = 1 — > Propagation 


C = 1 — > Propagation 


honeycomb 


C = 1 ^ Class IV 


C = 1 ^ Class IV 


quasi 


C = 1 ^ Class IV 
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